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RTCC REQUIREMENTS FOR MISSION G: SELECTING AND 

VERIFYING USBS DOPPLER DATA SOURCES DURING LM ASCENT AND DESCENT 

By Alan D. Wylie 


SUMMARY AND INTRODUCTION 

This note describes a program for computing MSFN high-speed Doppler 
tracking residuals used in the Real-Time Compute Complex (RTCC) for 
verifying and selecting USBS Doppler data sourt during LM ascent and 
descent. These same residuals are used to determine whether the primary 
guidance and navigation control system (PGNCS) or the abort guidance 
system (AGS) is the more correct LM navigation system if the two signifi- 
cantly disagree. Rendezvous radar observations are also processed to com- 
pute rendezvous radar residuals which provide an independent test on the 
accuracy of PGNCS versus AGS. The program has the capability of simul- 
taneously processing nondestruct Doppler data for up to four USBS ground 
stations tracking the LM during ascent and descent 9 one of which may be 
the transmitter. It is assumed that the Doppler biases for these stations 
have been computed. 

In order to check the validity of the ground tracking data, Doppler 
residuals are computed for each tracker using common telemetered LM state 
vectors from either PGNCS or AGS. The "residuals" should agree to within 
some expected tolerance. Should one station disagree, it can he replaced. 

Doppler residuals are actually computed using both telemetry vector 
sources so that the computer controller may display either set of residuals 
and so that both sets are available for the flight control vector source 
selection. For each vector source an average is taken of the valid station 
residuals in order to indicate which system is providing the better data. 
The rendezvous radar residuals are computed and displayed for the same 
purpose. 

The information in this note supersedes that in reference 1. 


GENERAL PROCEDURE FOR SELECTING AND VERIFYING DOPPLER DATA SOURCES 

Residual Computations 

The main purpose of this program is to verify the Doppler data of the 
ground stations tracking the LM during ascent and descent. The most 
efficient method is to calculate a Doppler residual - that is, to find for 
each station the difference between observed and computed Doppler values . 
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From each tracking station, time-tagged raw Doppler count observations 
arrive at the RTCC at a 10-per^-second rate. The observations pass through 
a MSFN preprocessor routine which makes simple edit checks and stores the 
valid data for the main program. Telemetered powered-flight LM state vectors 
arrive at an effective rate of one observation every 2 seconds from each 
onboard navigation system. Free-flight vectors predicted by the analytic 
ephemeris generator (AEG) and the lunar stay vectors generated using the 
selenography c radius (r) , latitude ((f)), and longitude (x) of the LM arrive 
at the same rate during the appropriate stages of the program. The telem- 
etered vector preprocessor and the predicted vector preprocessor make simple 
edit checks on the corresponding vectors and store the valid values for 
future processing. Rendezvous radar observations are telemetered down in 
the form of raw unprocessed data at a rate of one observation every 2 sec- 
onds. This data must also pass through a preprocessor routine which in- 
cludes simple validity tests on the raw data. 

The program starts the residual computations process by selecting a 
recent telemetry vector and iterating to find the transmitter- vehicle (up- 
link) range and vehicle-receiver (downlink) ranges. A linear interpolation 
scheme is used for each station to find the "pseudoactual" Doppler count 
which corresponds in time to the receiver time of the vector. 


The program recycles and obtains a second vector from the same source. 
A position vector is then computed at the time of the second position vec- 
tor and temporarily replaces the second position vector in the residual 
computations. In this manner the program avoids a large random error in 
the computed Doppler due to the large quantization error in telemetered 
position vectors. 


This computed vector is then processed to obtain uplink and downlink 

ranges and a value of Doppler count for each station. The observed and com- 
puted Doppler frequencies are calculated and subtracted to obtain the re- 
sidual. The residuals are converted from cycles per second to feet per 
second and plotted continuously in time as in figure 1. 



AY 3 
AY 2 

AYi 


Figure 1.- Station verification display. 
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Residuals will also be computed using the other vector source. Each 
set of valid station residuals is then averaged and plotted continuously 
in time on a flight control display to indicate which telemetry system is 
in better agreement with the MSFN trackers . 

Rendezvous radar range-rate residuals are also plotted on the flight 
control display. To compute one of the residuals, the program interpolates 
the current CSM ephemeris to obtain the CSM state vector at the time of the 
LM vector being processed. The CSM and LM vectors are then used to compute 
the range rate between the two vehicles . The smoothed value of rendezvous 
radar range rate evaluated at the time associated with the LM state vector 
is then used together with the computed range-rate value to construct the 
residual which is plotted on the display. 


Operational Modes and Procedures 

The ascent-descent program is designed primarily to function in the 
"average g" or powered- flight mode. However, it has been extended to 
cover short periods of free flight and the LM's stay on the moon. 

Program use actually begins approximately 5 minutes before ignition 
for powered descent and is initialized by MED 1. An RTCC predictor, the 
analytic ephemeris generator (AEG), will use the best LM anchor vector to 
project LM state vectors which are processed in the high-speed mode in the 
same manner as the powered- flight PGNCS and AGS vectors. Residuals are 
computed and displayed for each station for checkout purposes. When the 
onboard computer switches to the average g mode (approximately 30 seconds 
before ignition), the LM begins to send updated vectors. A switch is set 
at the RTCC console to initiate high-speed processing of actual telemetry 
vectors when they are received. . 

After touchdown on the lunar surface, the LM continues to send vectors 
used by this program for approximately 3 minutes. At the end of that time, 
a MED is entered to terminate the descent phase and thus cease residual 
generation. About 5 minutes prior to lunar lift-off, the MED to enter 
high-speed processing and MED 1 are reentered and the program resumes high- 
speed processing of RTCC-computed LM vectors and MSFN observations in prep- 
aration for the ascent phase. A switch will initiate the use of actual 
telemetry vectors when they are received. 

During the ascent stage, high-speed processing of PGNCS and AGS vectors 
continues until averaged vectors becpme available following engine cutoff. 

At that time the controller has the option either to continue using actual 
telemetry vectors or to switch and use the RTCC-predicted vectors based on 
the computed average vectors. If the average vectors are not available, 
the last powered- flight vectors will be considered in the option. The 
program continues to process the vectors (either actual or predicted) in 
the high-speed mode until an orbit switch is thrown, signifying the end 
of the ascent phase. 
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DETAILED EQUATIONS AND LOGIC 


Telemetry Vector Preprocessor 

Two telemetry vector preprocessors (appendix A) are used in this 
program in order that the telemetered powered- flight vectors and the pre- 
dicted free-f light or lunar-stay vectors may be treated separately. Each 
preprocessor is a totally independent routine which is started by executive 
control. Only one of the two may function at any given time. 

The telemetry vector preprocessor receives the LM powered- flight state 
vectors telemetered from the onboard PGNCS and AGS navigation systems in 
units of earth radii and earth radii per hour and performs simple tests on 
the vectors and time tags. Basic validity tests are included, (ref, 2), 
and all necessary granularities are applied. All vectors not edited out 
are stored in one table for future processing. The logic is described in 
flow chart A-l ( a ) . 

The predicted-vector preprocessor generates LM state vectors in the 
same units at the real-time rate while the LM is either in free flight or 
on the lunar surface. The logic is described in flow chart A-l(b). 


MSFN Preprocessor 

The MSFN preprocessor (appendix A) also is an autonomous routine 
started by executive control. It receives from the tracking stations high- 
speed data consisting of raw Doppler counts accumulated from count initiali- 
zation and not a count.', difference over a fixed time interval. Continuous 
tests are performed bo check if the "ascent-descent sites" MED has been 
entered. The change option of the MED causes one or more new trackers to 
be entered for which the appropriate storage areas must be reinitialized. 

The start option causes a new transmitter to be entered which results in 
all storage areas being reinitialized. Simple validity tests (ref. 2) and 
a test on the vector time tag are used to detect invalid data. Any nec- 
essary granularities are then applied. 

The preprocessor logic then proceeds to the edit routine. The initial 
two values of Doppler count for each station are stored, tagged invalid, 
and utilized in a linear extrapolation scheme to predict the Doppler count 
at the time of the third observation. The routine then tests to determine 
if the observation is greater than the last accepted value and less than 
some maximum value which is a function of the extrapolate . fcime interval. 

If the observation passes the test, it is stored but tagged invalid, and 
the program uses the last two stored values to linearly extrapolate to the 
time of the next observation. If the observation fails the test, however, 
it is discarded, and the routine uses the same two data points in repeating 


P 
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the extrapolation procedure. Once four successive data points pass the 
edit scheme, the four values and succeeding accepted values are stored 
oj^d tagged valid. Four successive failures, however, cause the routine 
to -e initialize. The program always tests if the linear extrapolation 
time interval is' less than 0 T< If not, the time interval is too long to 
linearly extrapolate, and the program must reinitialize. 


N 


The MSFN data passing the edit routine is stored into the arrays 
t . , where i represents the station and j represents the position 


j/i’ v j/i’ 

in the table. Flow chart A-2(a) describes the logic for the preprocessor. 
Flow chart A-2(b) describes the edit routine logic. 


Rendezvous Radar Preprocessor 

The rendezvous radar preprocessor (appendix A) is started by executive 
control. Input consists of raw unprocessed destruct Doppler count data and 
the associated time tag at the middle of the count interval. The routine 
performs elementary validity tests (ref. 2) and time tests before applying 
granularities (ref. 3) which convert the raw data to range rate in earth 
radii per hour. The routine then tests if the data is within a determined 
extreme interval. All observations passing the above tests are stored for 
the main program. The logic for this preprocessor is described in flow 
chart A- 3. 


Main Program Control Logic 

The main program control logic (appendix B) completes the list of 
processors which are started by executive control. It is broken down into 
five parts: the main program logic or driver, the iteration subroutine, 

the interpolation subroutine, the MSFN residual computation subroutine, 
and the rendezvous radar residual computation subroutine. Flow charts B-l 
through B-5 present the equations and logic for each of the five parts of 
the program. 

All variables must be initialized to zero at the start of the program. 
The orbit determination constants described in appendix B and the RNP ma- 
trix of current time are then read in. In addition, station characteristics 
and a lunar ephemeris must be available.' Tables of time and edited Doppler 
N-count observations are generated for n stations , 

The program begins each residual computation cycle by picking up a 
telemetry vector and waiting until the vector is T1 (time interval from 
vector receipt time to current time) seconds old in regard to current time. 
The first computations involve finding the ranges from the vehicle to the 
transmitter and receiver for each tracker. Since there is one transmitter 
and four receivers, one uplink range and four downlink ranges are computed. 
The iteration routine computes these ranges and the receiver- times for each 
station. 
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For each station the raw observations ax'e set up in the following 
arrays: t . ^ and ^ where i = 1, ...» n and j - 1, > • • , Ml. The in- 


dex Ml represents the most current observation at a station and is deter- 
mined by the time length of the data table. The time test is then performed 
for station i to determine if the receiver time obtained from the iteration 
routine lies within the time range of the current MSFN observation table. 
Failure of the test causes the program logic to check the next station. 

If all stations fail the test, the program logic transfers to the rendezvous 
radar residual routine. Otherwise it proceeds to the interpolation routine. 


The linear interpolation routine is then called to yield the inter- 
polated value of the Doppler count for each station at each receiver time. 
The main program then stores all necessary values derived from the onboard 
telemetry vector and proceeds to the rendezvous radar subroutine. 


Output from the rendezvous radar residual subroutine consists of a 
range-rate residual. It is computed vising a PCNCS or AGS telemetry vector, 
a free-flight CSM ephemeris , and the current table of rendezvous radar 
observations. The residual is plotted on a flight control display, and 
the program recycles to pick up the second telemetry vector necessary for 
the Doppler residual computations. 

When the second telemetry vector is receiveu from a vector source, 
tne program computes the position vector referenced to the same time. The 
program uses the position and velocity vectors of the previous telemetered 
LM state vector and the velocity vector of the current telemetered state 
vector (both vectors from the same source) in the computation and assumes 
linear motion of the vehicle. The position vector is computed in order to 
prevent a large quantization error from appearing in the Doppler MSFN re- 
sidual computations. The program uses the computed position vector in 
order to compute the second set of uplink, and downlink ranges (iteration 
routine) and values of Doppler count (interpolation routine) to be used 
in the residual computations. The logic then transfers to the MSFIT re- 
sidual routine. 


(This scheme for predicting a position vector using a previous 
position vector and the average of two velocity vectors is used only 
when actual telemetry vectors are being processed during the average 
g mode. During the non-average g stages, the nomine 1 vectors ar^ 
processed as they are received. ) 

The MSFN residual routine is provided with values of uplink and down- 
link range obtained from one telemetered vector and one computed vector 
from the same source and the two values of interpolated counts from each 
cf the n stations. The residuals are computed and displayed as shown in 
figure 1. The average residual 6Y computed using vector? from one vector 
source is then compared on a display with the average residual of the other 
vector source for the selected stations. The program then returns to pro- 
cess the actual telemetry vector vhc^e time tag is the same as the computed 
vector just processed. 
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Iteration Subroutine 

Input to the iteration subroutine (flow chart B~2) is the LM 
position vector, R(t^), in a selenocentric coordinate system with the 

X-Y plane parallel to the earth's mean equatorial plane, the X-axis through 

Aries at the beginning of the nearest Besselien year (KEY), and the Z-axis 

directed along the mean pole. Also input arc a vector time tag, t T „, ref- 

LM 

erenced to Greenwich midnight prior to launch, station characteristics for 
station i (i "1, n) , and a lunar epbemeris. The first computation 

transforms the mean NBY selenocentric vector to a mean NBY earth-centered 
inertial (ECl) vector by 


R'(t LM ) = R( 


V 


% 


where R„ is the position of the moon in the earth-r ferencod mean NBY 
M 

system. The next computation involves finding the range between the 
transmitting station and lunar orbiting vehicle. 

In computing the range, the program must iterate to find the time, 
t^, the signal was transmitted from the transmitting station in order 

to havi arrived at the vehicle at t^ (ref. 4)* Let c be the velocity of 

light, R^fc^) the geocentric inercial position vector of the vehicle at 

t T „ in the mean NBY coordinate system, and R (t*) the geocentric inertial 
LM s 

position vector of the transmitting station at the computed transmission 

time , t ’ . Then 


t' 



( 1 ) 


where 

P » iRMt^) - R s (f)| • 

For an initial gress of R q ( t f ) » we set t' = t^} thus, 

p - Ir’Uj^) - R^tjjPl . 


(> 2 ) 

(3) 


This p from equation (3) is then used in equation (l) for a new 
estimate of t*. The position of the station is then computed for the 
revised t' , a new p is calculated in equation (2), and a new t ' is calcu- 
lated in equation (l). This iterative process is repeated until the dif- 
ference between successive values of t' is less than or equal to 10- seconds 
or until six iterations have been made. 
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A refraction correction is then applied to the computed uplink range 
p to form p' in the following manner. (See ref. 5 .) Let 

r s = |R s (t')| 


and 


Then 


5 “ E,( V - E s (t,) 


p* = p + 


c , r p 
nk s 

R s (t') • p 


ik) 


where c^ is the product of refraction modulus and radio refractivity 

x 10 -6 , which are available from the station characteristic table. (See 

appendix B.) The uplink range, p', will then be stored as P™,. 

JLK 

The suoroutine then computes downlink ranges and times for the 
n stations successively. The only change from the uplink computations is 
in the time- iteration equation which now becomes 



(5) 


where t' is the computed time at which the signal should be received at 


the station if it left the vehicle at t, 


LM* 


After the refraction correction 


has been applied to a downlink range, p* is stored as P^/^ and t' stored 

ast (i=l, ...,n) and the subroutine picks up the proper station 
R/i 

characteristics for the downlink computation for the next station. 


The position of the station, R (t'), is computed for each successive 

s 

iteration of time, t' , from the station characteristics as follows: 

p 


' Let 0 be the angle between the true position of Greea- 

I wich at midnight prior to launch and Greenwich at t*. 

\ Then 

1 

0 = 0)t 1 

! 

i where o> is the rotational rate of the earth, 

j Then the angle, V, of the station is given by 

j X* m X + 0 (7) 
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where X is the longitude of the station. Finally, 

R (t ' ) is given hy 
s 


R (t ! ) a (RHP) 
s 


cos A f (r cos <f> f + h cos (f>) 
sin A*(r cos + h cos <f> ) 
r sin <J> ’ + h sin <f> 


( 8 ) 


All the above variables are defined in appendix B. 


Interpolation Subroutine 

The interpolation subroutine (flow chart B-3) is called to find a 
pseudoactual Doppler count, Vi- for each station at the receiver time 

t R ^ which is not likely to be a recorded observation time. Inputs to the 

subroutine are t R ^ (obtained from the iteration routine) and n tables of 

Ml observation times, t^^ , -versus Ml observed Doppler counts, N^, 

(i = 1, ..., n and j = i, ..., Ml). All stations for which the flag MSW(i) 
is zero are processed. 

The interpolation scheme tests the Vi of each station against the 

corresponding table of recorded observation times to find between which 
two recorded valid observation times Vi lies and between which two 

recorded valid observations lies. Once these two values of 

tj^ are found, the routine tests tj^ and t^ + ^^ to see if 

Vl/i - Vi - a T • 

If not, the time interval is too long to linearly interpolate. There- 
fore, do not interpolate, but set N R ^ equal to a dummy negative value 

denoted by a^, and start the time tests for the next station. If, however, 

Vl/i “ Vi ' 1 °T * 


the time interval is sufficiently short to accurately interpolate for 


N, 


R/i* 


Therefore, the following equation can be solved for N. 


’R/i* 


N, 


r/i ■ Vi + - ♦ 


>l/i u j/i 


m. (t _ t ) 

' n/i VP 


(9) 
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where i is the station number. The outputs N R /^ from the interpolation 
subroutine will then be passed to the main program. 


MSFN Residual Subrout ine 

Inputs to the MSFN residual computation subroutine (flow chart B-U) are 
two receiver times, two uplink ranges, two downlink ranges, two interpolated 
N-count observations, and a table of bias values. The subroutine will be 
called in two different places in the main program depending on whether 
the vector source is PGNCS or AGS; this selection is controlled by the LSW 
flag (see appendix B). The inputs to this subroutine can be controlled 
by the arguments of the calling statements (or equivalent); for example, 
for a PGNCS vector we might have: 

CALL MSFN RESIDUAL (Pt R ^, t R ^, Pp TR , Ppj^* P TR > P R /i> PN r/i> W R/i> 

AY , 6Y) 

Then for an AGS vector we would have: 

CALL MSFN RESIDUAL (At R ^, t R ^ , Ap TR , Ap R ^ , P TR , P R ^ j AN R ^ , N R ^, 
AY i , <5Y) 

Once the input arguments have been established, the program processes 
the residual subroutine, which is of the following form: 

SUBROUTINE MSFN RESIDUAL (t^., t^, p x# p 2/ ., p 3# p^, N^., N^, 
AY i , 6Y) 

Before computing a residual for an observation station, the validity 
of the two Doppler counts and N^ must be verified. A negative 

value for either N^ or N^ causes the residual computation of that 

station to be omitted, and the routine checks the Doppler counts for the 
next station. 


Let p 2 ^. and p^^ be the downlink ranges computed for the receiver 
times tjyj and t ^ and let p x and p 3 be the corresponding uplink ranges. 

The computed Doppler observation is obtained from the following equation 
(ref. h ). 


un.f, 


D t - (» 3 + V + 


Vtr 


2/i " *1/1 


[ 


( P o + Pjj/^) “ + Po/iH (W) 


2 /i 


>] 


\ 

} 




ft 
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where the constants co^, , c, b^ are explained in appendix B. A 
pseudoactual Doppler observation, S^, is obtained by 


S 


i 


N 2/i - Vi 
t 2/i " Vi 


( 11 ) 


and the residual is computed from 

AY. = (S i - D.) K, 


( 12 ) 


where K, the conversion factor from Hz to fps, is given by 


K = c. . 

2w 4 f TR 

where f is the frequency transmitted to the LM. These AY. are plotted 

TR 1 

as in figure 1. The average of the residuals from n' stations is then 
obtained by 


6Y = 


n' 


£ AY * 


&=i 


(13) 


where n' is the number of stations selected. This 6Y' will then be plotted 
against 6Y computed with the other vector source and displayed. Control 
is then transferred to the main program in order to process an actual 
telemetry vector and thus initiate the next residual computations. 


Rendezvous Radar Residual Subroutine 

The rendezvous radar residual subroutine (flow chart B-5) provides 
an additional check for testing PGNCS versus AGS. Inputs to this routine 

are a table of range-rate observations versus time (p^, t^; RR = 1, M2), 

a table of selenocehtric CSM ephemeris state vectors versus time (R„ om , 

V CSM» t CSM^» ^<1 a telemetered state vector (Rj^, Vj^, t^). A routine 

is called to obtain specific rendezvous radar observations in the time range 
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from t -20 seconds to t T „ + 6 seconds (see appendix C). If less than 
LM LM 

three such observations exist, the rendezvous radar computations are omitted 
for the particular telemetry vector. 

If sufficient rendezvous radar data is available, a sixth-order 
Lagrange interpolation routine interpolates the CSM ephemeris table to 
obtain the CSM state vector at the time of the LM state vector being pro- 
cessed. The rendezvous radar range-rate data is smoothed by a second- 
degree, least-squares filter (see appendix C). This second-degree curve 
is then evaluated at the time of the telemetry vector to obtain the de- 
sired range rate observation. 

Tbs position vectors of the LM and CSM allow the program to compute 
the range between the two vehicles using the equation 


p C0MP = R CSM " R LM' 

Then the range rate can be computed by 

_ P C0MP , 

p C0MP .- | * V CSM ' 111 

|p C0MP' 


(lit) 


(15) 


This computed range rate is then compared to the observed rendezvous 
radar range rate, and the residual, Ap, is displayed on the flight control 
display in units of feet per second. The logic then, returns to the main 
program to pick up a new telemetry vector. 


CONTROLS AND DISPLAYS 


The following are general descriptions of the various controls and 
displays used in this program. Reference 2 presents detailed; descriptions. 

Manual Entry Devices 

1. Ascent-descent sites .- This MED specifies up to fpur.- sites to he 
processed in the high-speed mode. This MED also allows the controller to 
change any of the sites being processed As well as to reinitialize if a 
change in the identity of the two-way site occurs. 
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2. Ascent-descent scaling factor .- This MED allows the controller 
to adjust the ordinate scale for either the MSFN or rendezvous radar resid- 
ual displays . 


Push Button Indicators 

1. Sources to be considered in the average .- This group of five 
PBI's determines which sites are to be used in the MSFN average residual 
computation. The ability exists to edit sites in or out of the average 
computation. 

2. TM selected source .- This group of two PBI's determines whether 
to display the MSFN residuals computed using the PGNCS vectors or the MSFN 
residuals computed using the AGS vectors. 

3. Processing mode .- This group of two PBI's determines whether 
actual telemetered vectors or AEG-predicted vectors are to be used in the 
MSFN residual computation. 


Displays 

1. MSFN source validation 1 and 2 .- This display provides an analog 
trend curve of range rate versus time for the residuals of two MSFN sources 
computed using PGNCS or AGS telemetered vectors. 

2. MSFN source validation 3 and 4 .- This display provides the same 
information as display no. 1 for two other MSFN sources. 

3. TM source comparison .- This display provides an analog trend curve 
of range rate versus time for both telemetry sources. Average MSFN resid- 
uals and rendezvous radar range-.rate residuals will both be displayed. 


PRECEDING PAGE 


blank not filmed, 
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APPENDIX A 

FLOW CHARTS FOR PREPROCESSOR ROUTINES 
TELEMETRY, MSFN, AND 


RENDEZVOUS RADAR DATA 
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Variables Used in LM State Vector Preprocessors 


t 

a 


time tag (g.e.t.) of a PGNCS or AGS telemetry vector 
unprocessed by the preprocessor routine. 


|R*| 

|V*| 

t*. 


selenocentric position and velocity of a telemetered 

vector at t in mean NBY coordinates 
a 

maximum acceptable value for | R j » 80 earth radii 

9* 

maximum acceptable value for |V | , 50 000 fps 

9- 

current time plus 6 seconds 


^MIN 


current time minus 20 seconds 


t LM 

V( v 

LSW 


. time tag of a current PGNCS or AGS telemetry vector ref- 
erenced to Greenwich midnight prior to launch 

selenocentric state vector (PGNCS or AGS) of the vehicle 

at time t T „ in mean NBY coordinates 
LM 

switch set by the identification code of the telemetry 
vector: LSW » 0 for a PGNCS vector, LSW 4 0 for an AGS 

' vector , LStf = 0 for all nominal vectors. 
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Flow chart A-l,- Preprocessor for LM vectors. 

(a) Logic for the telemetered vectors - Concluded. 
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Cl(i) 


C2(i) 


C3(i) 


W k/i 


Vi 


MIN 

t* 


V/i 


V i 


T 

NMAX 


HIGH 


Variables Used in MSFN Preprocessor 

station counter of Doppler data which do not pass the 
extrapolation test; when Cl(i) equals four, all counters 
for station i are reinitialized and the extrapolation 
edit routine starts over 

station counter of Doppler data which pass the extreme 
count interval test; C2(i) = 3 initializes the use of 
extrapolating for Doppler 

station counter of Doppler data which pass the extrapola- 
tion test: C3(i) >_ k causes the data to be tagged valid 
when stored by the preprocessor routine 

an array of raw, non-destruct Doppler counts recorded at 
station i; the element N^y ^ represents the current value 


time array for the N^y^ 
current time minus 20 seconds 


current time plus 6 seconds 


an array of raw, non-destruct Doppler counts 
stored by the MSFN preprocessor routine; the 

represents the last accepted value of N^ 


accepted and 
element N^ 


time array for the N 


<j/i 


maximum time interval allowed for linear extrapolation 
of Doppler 


maximum acceptable count difference between the actual and 
extrapolated values of ; UMAX ts a function of the time 
interval t J/;1 - t k/± 


maximum acceptable value of 
t *e requirements have been 


Vi 

met 


tested when all previous 




Flow chart A-2.- MSFN preprocessor, 
(a) Logic for the preprocessor. 
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(a) Logic for the preprocessor - Concluded. 





Flow chart A-2.- MSFN preprocessor, 
(b) Logic for the edit routine. 
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Variables Used in Rendezvous Radar Preprocessor 


time tag of a rendezvous radar observation unprocessed 
by this routine 

raw unprocessed rendezvous radar destruct Doppler count 
which represents the count of a frequency comprising both 
the Doppler frequency and a bias frequency over a time 
interval 

current time minus 20 reconds 
current time plus 6 seconds 

scale factor required to obtain the range rate in earth 
radii per hour. It is of such polarity as to make 

• 

positive for increasing range, 
range-rate bias frequency 

counting interval over which is observed 

time tag of a rendezvous radar observation processed 
by this routine 

lower limit for acceptable rendezvous radar range-rate 
data 

upper limit for acceptable rendezvous radar range-rate 
data 

processed rendezvous radar range-rate observation at t^ 





m 
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Flow chart A-3.- Logic for the rendezvous radar preprocessor - Concluded 
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FLOW CHARTS FOR PROGRAM CONTROL LOGIC 
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Variables Used In The Main Program 




v( V 


selenocentric LM state vector (PGNCS or AGS) of the vehicle at 

time t T „. in MNBY coordinates 
LM 


R(tj^j) LM selenocentric position vector (PGNCS or AGS) used in order 

to lessen a quantization error in the Doppler residual computa- 
tions 

PSW switch used to identify the PGNCS vectors; PSW 4 0 denotes an 

actual telemetered PGNCS vector, PSW = 0 denotes a computed 
PGNCS vector 

ASW switch used to identify the .AGS vectors ; ASW 4 0 denotes an 

actual telemetered AGS vector, ASW = 0 denotes a computed AGS 
vector 

T1 elapsed time from vector receipt time to current time 

T2 time increment used in delaying vector processing 

C4 counter which counts the number of stations that cannot meet 

the time tests for a particular telemetry vector 

MSW(i) switch which denotes MSFN data: MSW(i) 4 0 implies time invalid 

data for station i; MSW(i) = 0 implies valid data which can be 
processed in interpolation routine 

j position number in an observation table 

Ml number of MSPN observations copied into a data table 

KSWP switch which denotes the PGNCS vector: KSWP = 0 for the first 

PGNCS vector, KSWP 4 0 thereafter 

TMAX the maximum time interval between two successive telemetry 

vectors from the same source for which a Doppler residual may 
be computed 

KSWA switch which denotes the AGS vectors KSWA = 0 for the first 

AGS vector, KSWA 4 0 thereafter 
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i 


PI>( V' 

pv(t M ) 

^LM 

n R/i > 
Pp R/i 
Pp TR 
PN R/i ■> 


temporary storage of previous values of t LM’ 

tR/i’ p R/i’ P TR’ and N R/i’ res P ectivel y » resulting from a. PGNOS 
telemetry vector 


1 



temporary storage of previous values of R(t^^, V(tj^) , t^ 
tR/i’ p R/i » P TR » and N r/i’ respectively, restating from an AGS 
telemetry vector 



GR(t M ) \ generalized temporary storage locations for Pt^, PRCt^), 
GV(t ) I PV(t LM ; or At M’ AY(t Uf) in the com P utation of 6( V 
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Flow chart B-l.- Logic for the main program 
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OBSERVE TABLE 
OF LAST Ml MSFN 
OBSERVATIONS 


IS t R/ . WITHIN TIME 

1 1— 

SPAN OF CURRENT 
MSFN TABLE OF 
STATION i? 




<e 

^ t l/i St R/i St Ml/i 




i = i + 1 
MSW (i) = 0 


MSW (i) =1 
C4 = C4 + 1 


ALL STATIONS 
PROCESSED 


ARE THERE ANY MSFN 
STATIONS WITH VALID 
DATA WHICH CAN BE 
PROCESSED BY 
INTERPOLATION SCHEME? 


ASW= OH — <ASW= Q> 
IYES 


:lsw=o. 


ASW = 1 
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1uominal> ^ yes 

s VECTOR/ / ^ 

Jno 

[psw=ti 


Flow chart 3-1.- Logic for ths main program - Continued. 
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Variables Used in the Iteration Subroutine 


R,( V 


i 

ISW 


LOOP 

P,P ' »t' , 
+" 


0 


0 ) 


ECI position vector of the vehicle at t^ in the Aries mean 
NBY coordinate system 

position vector of moon, ECI coordinates 
station number, where i = 1, n 

switch which denotes uplink or downlink range computation: 
ISW = 0 for uplink range computation, ISW $ 0 for downlink 
computations 

counter which limits the number of iterations 
temporary storage locations used in the iteration cycle 

angle between true-of-date position of Greenwich and the 
position of Greenwich at midnight prior to launch 

rotational rate of the earth 


DT difference between UT1 and UTC , DT = UT1 - UTC 

A f A + 0 

A longitude of station 

R (t 1 ) ECI position vector of a station at the computed time t f 

3 in the Aries mean NBY coordinate system 

R (t ) initial guess of R (t 1 ) 

S JjM s 

T 

(RNP) matrix which converts from the Greenwich position at mid- 

night prior to launch true-of-date coordinate system to 
the Aries mean NBY system 

The following constants are found on the station characteristics 

tape: 


r cos <j>* + h cos <(> 


projection of position vector of the station 
into the x-y plane in the ECI coordinate system 


r sin <t>* + h sin <)> 


z coordinate of station in ECI system 
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P 



t TR 

P TR 

Vi 

P R/i 


where r radius to ellipsoid below station 

(f> ' geocentric latitude of station 

4) geodetic latitude of station. 

h height of station above ellipsoid 

vector difference between R’(t T ,.) and R (t') 

LM s 

magnitude of the vector R (t') 

s 

product of the radio refract ivity x 10~ 6 and the decay 
constant calculated from radio .refract ivity with an ex- 
ponential reference atmosphere; both variables are found 
on the station characteristics tap? 

time the signal must be transmitted from the transmitting 

station to arrive at the vehicle at t T „ 

LM 

uplink range from the transmitting station at t' to the 

vehicle at t T „, 

LM 

time at which the signal would arrive at station i if it 

had left the vehicle at t T „, 

LM 

downlink range from the vehicle at t T „„ to station i at 
. / . i„ \ LM 

Vi (l = 11 nl 

number of stations 




I 

1 


1 

j 

i 

\ 


Flow chart B-2.- Logic for the iteration subroutine, 











43 


Variables Used in the Interpolation Subroutine 


T 


Vi 


dummy negative value for Vi which signifies that a 

residual will not he computed for station i for the current 
telemetry vector being processed 

the maximum time interval over which the Doppler data may 
he linearly interpolated 


J/i 


and N 


an interpolated N-count situated between N 


0+1/ i 







N l/i 

N 2/i 

D. 

1 

“3 


“4 


b. 

i 


f 


TR 


c 



P 3 I 

p 2/i 

P 4/i 

S i 

AY. 

l 

K 

STA(i) 

(SY 

n' 


45 


Variables Used in the MSFN Residual Subroutine 


interpolated values of N-count at t^^ and t^^ 
computed Doppler frequency at station i 

a bias introduced by the tracking equipment that guarantees 
that the Doppler shift will never become negative 

a multiplying factor which adjusts the frequency of the 
signal at the spacecraft 

bias for station i on the actual Doppler observables 
(previously solved for in a preascent routine) 

frequency transmitted to the LM» as specified on the 
station characteristics tape 

velocity of light 


successive computed observation times for station i 


uplink ranges computed for and t 


respectively 


downlink ranges computed for t^^ and ^2/i ^spec^ively 


observed Doppler frequency shift over tg^ - 
residual computed for station i (actual minus computed) 
constant which converts Hz to fps 

flag subject to external control via the station PBI 


average value of the AY^ 

number of stations used in computing 6Y; the value of 
n' is subject to external control via the station PBI 



INPUT: ORBIT DETERMINATION CONSTANTS; 
TABLE OF DOPPLER BIAS VALUES; THE * 
FOLLOWING VARIABLES OBTAINED FROM 
THE ITERATION AND INTERPOLATION 
SUBROUTINES: t 2/| , P y P 2/[l 

Py ^ 2 / 1 ' ^ ~^’ i 


"START" 


DO NOT USE DATA 
FROM STATION i 
IN COMPUTATION 
OF 8Y 


0<N l/i <N 2/i 


COMPUTED "DOPPLER FREQUENCY" 

TR l 

n — cu t h ! - K 

[ ( P 3 + /> 4/i > " ( ^1 + P 2/& 

U 1 3 + D | + eft -t ) 

1 i i Clt 2/| tjyjJI 


ACTUAL "DOPPLER FREQUENCY" 

_ N 2 /r N l/i 

r t — 

”2/1 r l/I 


| COMPUTE RESIDUALS j 

li 

< 

M 



DISPLAY 

AY, 


STA (i) 
=0 


Flow chart B-4.- Logic for the MSFN residual subroutine. 
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Variables Used in Rendezvous Radar Residual Subroutine 


p s' 


time tag of first rendezvous radar observation after 

rendezvous radar range-rate observation used in smoothing 
routine 


time tag of p r 
s 


CSM 


V 


CSM 


CSM ephemeris state vector whose time tag is the time of the 
LM telemetry vector being processed 


LM 

/ 

m number of rendezvous radar observations used in the data 

smoothing subroutine 

h counter of rendezvous radar observations 

At, time interval between the time tag of a rendezvous radar 

1 value in the table and the time tag of the telemetry vector 
being processed 

p, raw rendezvous radar range-rate value which corresponds to 

h . , 

At,_ 


\ rendezvous radar data summations for m points 


7 . 

MAT 1 
MAT 2 


matrices which define the equations used in the routine 


MAT 3 



^9 

coefficients of the second-degree curve which best fits the 
m rendezvous radar values being processed 

# 

P^-do rendezvous radar range-rate value observed at t TM 
UJdo IjM 

°COMP computed range from LM to CSM 

P COMP magnitude of p COMP 

p COMP computed range rate from LM to CSM 

Ap range-rate residual (actual minus computed) 
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Flow chart B-5,- Rendezvous radar residual subroutine, 

(a) Logic for the subroutine - Concluded* 









m 



Row chart B-5,- Rendezvous radar residual subroutine. 

(b) Logic for obtaining data for data smoothing routine - Concluded. 




INPUT: TABLE OF 
m RENDEZVOUS RADAR 
OBSERVATIONS, TIME TAG 
OF LM TELEMETRY VECTOR 
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Flow chart B-5.- Rendezvous radar residual subroutine 
(c) Logic for the data smoothing subroutine. 
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Flow chart B-5.- Rendezvous radar residual subroutine, 
(c) Logic for the data smoothing subroutine - Concludes 
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PROCEDURE FOR SMOOTHING RENDEZVOUS RADAR DATA 


Data Acquisition Routine 

Rendezvous radar data is received via the telemetry downlink in the 
raw unprocessed form at a rate of one observation every 2 seconds* It is 
initially smoothed by a series of validity tests as described in reference 2. 
Further smoothing is necessary, however, as is an interpolation scheme for 
the data. A second-degree, least-squares curve fit has been selected to 
accomplish both purposes. 

The objective of this routine is to find six rendezvous radar range- 
rate values in the data table which are in the time vicinity of the time 
tag of the telemetry vector being processed (t^) . This time interval is 

from t T .. - 20 seconds to t TW + 6 seconds. The rendezvous radar residual 
LM LM 

routine can function if only three, four, or five such points are available; 
otherwise computations are omitted for the particular telemetry vector. 

The logic for this routine begins by locating the two data values 
between which the time tag of the telemetry vector lies. The program then 
goes forward in time from t^ to t^ + 6 seconds or to the end of the table 

and records all data values (not more than 3) in the interval. It then 
goes back in time from t^ and records all data values until a total of 

six values has been recorded or until it reaches the time t T *. - 20 seconds 

LM 

or the end of the table. 

For each data value the time interval between its time tag and t^ 

is also recorded. The logic then returns to the rendezvous radar residual 
routine. 


Data Smoothing Routine 

Any second-degree curve of range rate versus time can be represented 
by the equation 


p = + a^At + a^At 2 . (l) 

Let the x’s below represent the raw data points p ^ 9 •••» P ffi 

acquired from the data acquisition routine to which it is desired to fit 
the second-degree least squares curve. 


a 


a 



t 


For each data point the difference between the raw observed value and 
the corresponding value on the curve is its residual. The objective of 
this routine is to solve for the coefficients a,Q, a^, a^ which determine 

a unique second-degree curve such that the sum of the squares of the resid- 
uals of the data points is a minimum. 

The zero time to which the At of equation (l) is referenced is the 
time tag of the telemetry vector being processed. Let t^ be this refer- 
ence time and let h be the counter of the data points recorded in the data 
acquisition routine. Then equation (l) can be evaluated at different 

points (t ,) on the curve as follows: 
s 


Pic " a 0 + a l At l + a 2 At l 
p 2c = a 0 + a l At 2 + a 2 At 2 

( 2 ) 


p mc " a 0 + a l At m + a 2 At A 



6 1 


where 


At, = t TVI - t , 
h LM s' 


p h = p s t for h = 1, . , m. 


The siM of the squares of the residuals is given by 

3 = j(aQ + a^At^ + a 2 At 2 ) - pj 2 + £(a 0 + a-jAtg + & 2 At 2 ) 

+ . . . + |(a_ + a.. At + a„At 2 ) - pi 2 
[ 0 1 m 2 m ' mj 

In order for S to be a minimum, 

3S 




( 3 ) 


9 a „ 


9 S 
9 a. 


= 0 


= 0 


W 


9 S 


9 a. 


= 0 


Solving the partials, we get three equations with three unknowns, 
which in matrix form can be expressed as 


[ a 0 a l a 2 ] 


It 

L h=1 


m 

m 




E 4t h 


m 

h=l 

h-l 



m 

m 

m 



E 4t » 

h=l 

h=l 

h= 

i 

m 

m 

m 


E 4t h 


h*l 

h=l 

h=l 

«■ 

m 

m 

- 



h p h £ 

p h 


h=l 

h=l 


• 


( 5 ) 




I 1 
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By applying a matrix inversion scheme to the 3 x 3 matrix, the co- 
efficients a^, a^, a^ can he readily computed. These coefficients deter- 
mine a unique second-degree, least-squares curve which can then he eval- 
uated for any At in the defined interval. In particular, at the time tag 
of the vector. 

At = 0 

» 

and the desired smoothed rendezvous radar range rate, p , is given by 

OBS 

P OBS = a 0 
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